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Abstract 

We present conjectured exact expressions for two types of correlations 
in the dense 0(n = 1) loop model onLxoo square lattices with periodic 
boundary conditions. These are the probability that a point is surrounded 
by m loops and the probability that k consecutive points on a row are 
on the same or on different loops. The dense 0(n = 1) loop model is 
equivalent to the bond percolation model at the critical point. The former 
probability can be interpreted in terms of the bond percolation problem 
as giving the probability that a vertex is on a cluster that is surrounded by 
[m/2\ clusters and \_(m + 1) /2J dual clusters. The conjectured expression 
for this probability involves a binomial determinant that is known to give 
weighted enumerations of cyclically symmetric plane partitions and also 
of certain types of families of nonintersecting lattice paths. By applying 
Coulomb gas methods to the dense 0(n = 1) loop model, we obtain new 
conjectures for the asymptotics of this binomial determinant. 

1 Introduction 

In this article we consider the dense 0(1) loop model [Hj defined on an L x oo 
square lattice with periodic boundary conditions. The geometry of the lattice 
is thus a cylinder with circumference L. The dense 0(n) loop model can be 
defined as follows. At each vertex the four edges meeting there are connected 
with equal probability in either of the two ways shown in Fig. ^ 

If L is even, the set of edges connected to each other form closed loops 
with probability 1. If L is odd there is a loop spanning the entire length of 
the cylinder. To a state consisting of I closed paths (loops) a weight of n l is 
assigned. In this article we consider exclusively the case n = 1, in which all 
configurations have equal weight. At this point the model can be mapped to 
the bond percolation problem at criticality, the six vertex model and the XXZ- 
spin chain at A = — ~ (see [B]). The corresponding bond percolation problem 
is defined on one of the sublattices of the dual lattice. The states of the dense 
0(1) loop model are in direct bijection with bond configurations of the bond 
percolation problem, see Fig. |2] 
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Figure 1: The two vertices of the dense 0(1) loop model 
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Figure 2: The mapping of a loop configuration to a bond configuration of the 
corresponding bond percolation problem. A bond is cither put on an edge of 
the square lattice formed by the • or on the dual edge orthogonal to it on the 
square lattice formed by the o. 

Recently, a number of conjectures have been obtained about correlations in 
this model. These involve probabilities that loops intersect a horizontal cut 
between vertices (henceforth referred to as a row) in certain prescribed ways. A 
row intersects loops at a total of L points. We define the connectivity state of 
the row as the way these L points are connected to each other by the loops via 
the half cylinder below the row. See Fig. [3| for an example. Connectivity states 
can be denoted by strings of parentheses |18j . If a point at position i is connected 
to a point at position j, then this is represented by a parenthesis at position 
i matching with a parenthesis at position j. Because of periodic boundary 
conditions, a right parenthesis at position i can match with a left parenthesis at 
position j also if i < j. For odd L, one point will not be connected to any other 
points on the row. That point is on an unpaired line. This can be denoted by 
a line | . For example, the connectivity state ) | ( (), denotes the state in which 
the first point is connected to third point over the boundary, the second point 
is not connected to any other point, and the fourth point is connected to the 
fifth point. 

Several conjectures about the properties of this probability distribution have 
been obtained, see |7l 1221 |2"51 1121 118| . These results can be interpreted as conjec- 
tures for probabilities of certain events on the top row of a half infinite cylinder. 
A recent summary of these results is given in 12 . Surprising connections with 
the enumeration of certain symmetry classes of alternating sign matrices |24] 
have been discovered. In this article we will make use of the following observa- 
tions: If the sum of all probabilities is normalized to the number of L x L half 
turn symmetric alternating sign matrices, denoted as Aht (L), then all proba- 
bilities are conjectured to be integers. It is conjectured that the probability of 
the connectivity state (((. . .))) for even L and | (((. . .))) for odd L with all the 
parenthesis nested inside each other is 1 in this normalization. For even L, it 



2 



V J V V_ 













J 






r ^ 
L J 




r 


( 

J 


r 








r ~\ 




f 


r > 






\ J 


j 











V J 



r 



v j 



V / 



r 



r 



v J 



( > 



Figure 3: Part of a typical configuration of the 0(1) loop model on a 6 x oo 
cylinder, periodic in the horizontal direction and extending to infinity in both 
vertical directions. The connectivity state of the half cylinder below the row 
indicated by the dashed line is () () (), above the row it is ) () ( (). 



has been proven JH] that: 

L/2-1 

9 i r 

4 (2fc + l)! 2 (2fc- l)! 2 



o L ^ r 1 3 (3fc + 2)!(3fc-l)!fc!(fc-l)! 
A HT (£)=2 [I —————— 2 (1) 



fe=i 

For odd L the number of L x L half turn symmetric alternating sign matrices 
is conjectured to be |24| : 

^ 4 (3fc)! 2 fc! 2 
A HT (i)- n (2) 

In some expressions for correlations in this paper the total number ofiixn 
alternating sign matrices, A(n) appears. This is given by |27l I15| : 

k=0 y ' 

A convenient way to compute probabilities of events at a row on the com- 
plete cylinder is to consider separately the connectivity states on both half 
spaces below and above the row. Since events above and below the row are 
independent, the probability distribution over the connectivity states contains 
all the information needed. 
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Tabic 1: Values for Q (L, m) 



m 


L 
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12 
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13167 


20048886 


247358122583 
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13446619579882992 
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Table 2: Values for Ci/ 2 _ m (L) 



m 


L 




2 


4 


6 


8 


10 


12 





3 


72 


13869 


21537270 


270190791369 


27414197906689626 


1 


1 


29 


6084 


9952274 


129464229047 


13499435968309125 
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351 


744193 


11416400590 


1338642053600985 
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112088583 


17605464402686 
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18952951005 
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2 Probability of being surrounded by m loops 

In |19| we have presented a conjectured exact formula for the probability P (£, m) 
that a face of the lattice on a cylinder of circumference L is surrounded by m 
loops. This conjecture is only valid for even L. Here we will explain how we 
obtained this conjecture and also generalize this conjecture to odd L. 

Since it is conjectured that all probabilities of connectivity states are in- 
tegers when the sum of all probabilities is normalized to be ^4ht [L), multi- 
plying P (L, m) by Aht (L) should be an integer. Let's define Q (L, m) = 
P (L, m) Aht (L) 2 - The values for Q (L, m) that we used to guess this function 
are given in tabled 

Guessing integer sequences is not difficult when the integers factorize into 
small primes, see e.g. In such cases simple product formulae can 

often be found by inspecting the prime factorizations of the integers. In this case, 
however, some integers have large prime factors. By inspecting the numbers 
Q (L, to), we noted that they are close to certain coefficients of the characteristic 
polynomial of the Pascal matrix. We denote by C p (L) the absolute values of 
the coefficients of the characteristic polynomial of the L x L Pascal matrix: 



det 

Kr,s<L 



r + s-2 
r-l 



— xS r 



= '£C n (L)(-x) n (4) 



n=0 

It turns out that Q(L,m) is approximately Cj,/ 2 -m (L) for m < L/2 — 1. If 
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to > L/2 — 1 then the two quantities are exactly the same: 

Q(L,L/2) = C (L) = 1 
Q(L,L/2-l) = C 1 (L) U 

In table|3we have listed the values of C^/ 2 -m (L)- By comparing the entries of 
tables ^and [21 we see that 

Q(L,L/2-2) = C 2 {L)-L/2 (6) 

It was a little more difficult to obtain this guess for Q (L, L/2 — 3): 

Q (L, L/2 - 3) = C 3 (L) - (L/2 - 1) C x (L) (7) 

The results obtained so far led us to the idea that Q(L/2 — k) might be a 
linear combination of Ck-2r (L) for r = 0... L^/2J, where the coefficients are 
polynomials in L. Using this idea we obtained: 

Q (L, L/2 - 4) = C 4 (L) - [L/2 - 2) C 2 (L) + | (j - 3 

Q(L,I/2-5)=C 5 (L)-(I/2-3)C 3 (I) + V^-4j - lj C x (L) 

(8) 

At this point we were starting to run out of data. The expression for 
Q (L, L/2 - 5) had been guessed using only Q (10, 0) and Q (12, 1). Only Q (12, 0) 
was left to guess Q (L, L/2 — 6). When we substitute L/2 — > to + k in the poly- 
nomials multiplying C r (L) in the expressions for Q (L, L/2 — fc), we see that 
the polynomials depend only on m. The results obtained so far led us to the 
conjecture: 

L/A-m/2 

Q(L,m)= K r ( m )C L/2 - m -2r(L) (9) 



where 



K (to) =1 

K x (to) = - (to + 2) 

K 2 (to) =i (to + 1) (to + 4) 



(10) 



Using the value for Q (12, 0), we can see that K% (0) should be —2. Since 
Ki(0) = -2 and K 2 (0) = 2, we guessed that isT r (0) = 2(-l) r for r > 0. 
Assuming that this pattern holds and that K r (to) factors into linear terms, we 
guessed that the sequence of polynomials could be — (m + 2), 1/2! (to + 1) (m + 4), 
-1/3! (to + 1) (to + 2) (m + 6), 1/4! (to + 1) (to + 2) (to + 3) (to + 8), etc. We 
thus put our faith in the formula: 

L/A-ra/2 , . 

Q(i,TO) = c i/2 _ m (L)+ V (-iyc L/2 _ m _ 2r {L)^--^ m + r (11) 

z — ' to + r \ r I 

r=l x ' 
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Besides fitting the data, any guess for Q (L,m) has to be consistent with the 
normalization: 

L/2 

Y,Q(L,m) = Am(L) 2 (12) 

TCI— 

Also, we had noted that putting a factor (—1)™ in this summation yields the 
total number of L x L alternating sign matrices: 

L/2 

Y j {-l) m Q{L,m)=A{L) (13) 

m=0 

We have proved that these two relations follow from iJUJ. See the next section 
for more examples of such identities. Also we have generated data for L = 14 
for Q (L, to) and verified that the formula predicts the correct values. 

For odd L it turns out that Q (L, to) is close to C(i_i)/2— m (L). In this case 
we were led to this conjecture: 

(£-1) m 

Q{L,m)= ^Vir 

r=0 

(14) 

The probabilities P (L,m) can be interpreted in terms of the related bond 
percolation problem, see Fig. Loops are boundaries of percolation clusters, 
separating clusters from dual clusters. If a point is surrounded by no loops then 
that point is on a cluster that wraps round the cylinder. A point surrounded 
by k loops is on a cluster that is surrounded by [(k + 1) /2J dual clusters and 
[k/2\ clusters. 



C(L-i) „ (L) — G th-i) 

- — ^ — L —m — 2r v ' - — o — -- 



■m — 2r- 



_x(£) 



3 Arbitrary weights for loops winding round a 
point 

Using the conjectured expression for P (L, to), we can give arbitrary weights to 
loops winding round a point. By (f>(L,z,x) we denote the operator that gives 
loops that wind round the point x a weight of z. By translational invariance its 
expectation value is independent of x, and in the following we will suppress the 
argument x. We can write: 



L/2 

>{L,z))= Y^P{L,m), 

m—0 



(15) 



Conjecture (|11|) and l|14|) imply that: 



>(L,z — 2 cos (9))) = — — ®\ if L \-, even. 



Aht (i) 



<(L,z= 2cos(6»))) 



D(L,i 



(16) 



2 cos (0/2) Ak T (L) e 



Here 



D(L,6) = exp(-i6L/2) det 

Kr,s<L 



r + s-2 
r-l 



if L is odd 



exp (id) <5 r . ; 



(17) 
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Table 3: Values for D (8) for 8 a multiple of 7r/3. Except for the first entry of 
the Table, the expression in terms of A (L) and Aht (L) is new. 
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even L 


odd L 





A HT (2L) 
ML) 


A HT (2L) 
ML) 


tt/3 




V3A HT (LY 


2tt/3 


A(L) 


A(L) 


TT 








The proofs of these relations can be obtained by substituting the conjectures 
for P (L, m) in l|15f) and performing the binomial summations. We show here 
how this works for even L. For odd L the proof is similar. For even L, 
implies: 

L/2 L/2 

Y,Q(L,m)z m = Y,B k (z)C L/2 _ k (L) (18) 

m=0 fc=0 

where 

So (z) = 1 

B fe (,) = g(-ir-^( fe - r )^for fc >l (19) 

The summand in this equation, which we name f(k,r,z), satisfies the recur- 
rence: 

/ (k + 2, r + 1, z) - zf (k + 1, r + 1, z) + f (k, r,z) = (20) 
Summing this recurrence over r from to k/2 yields: 

B k+2 (z) - zB k+1 (z) + B k (z) = (21) 

From (fmf> we see that B 1 (z) = z and B 2 (z) = z 2 — 2. The solution of J2IJ 
satisfying these two initial conditions is 

B k (z) = 2 cos (k6) (22) 

where 2cos(0) = z. Inserting this in l|18l) and using the relation C^^+k = 
C L /2-k El yields (H| for even L. 

Only for 8 a multiple of tt/3 are exact evaluations of D (L, 8) known |l()l IHl 
114) . For 8 = it is known that the determinant is given as Aht (2£) /A (L). We 
have found that for 8 any multiple of tt/3, D (L,9) can be expressed in terms 
of A (L) and Aht (L). These expressions are given in Table |3| To prove these 
formulas, one can divide the product expressions for these determinants given in 
|l()U14j by the expressions given in Table substituting the formulas for A (L) 
and Aht (L) given by Q, J5J and 10. The resulting expressions then become 
rational functions which are easily proven to be identically 1 . 
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4 Asymptotic behavior of (0 (L, z)) 



In this section we will obtain the asymptotic behavior of (4> (L, z)) by apply- 
ing Coulomb gas techniques 1201 IS] to the dense 0(1) loop model. First one 
randomly assigns an orientation to the loops of a state. One then assigns a Boltz- 
mann weight to the oriented loops so that the sum of the Boltzmann weights 
of the two orientations for each loop is unity. This is achieved by assigning 
a weight of exp(i7r/12) each time a loop takes a right turn and exp(— i7r/12) 
each time a left turn is taken. To a configuration of oriented loops a height 
function h (r) can be defined on the faces of the lattice as follows. Rotate all 
the arrows at the boundaries between the faces 90 degrees anti-clockwise. The 
rotated arrows define the (discrete) gradient of the height function. The differ- 
ence between two neighboring faces is ±1. It is assumed that in the scaling limit 
the discrete height function may be replaced by a real valued function with a 
Gaussian measure M. [{h (r)}] ~ exp (— H g (h)), where 

H E = ^Jd 2 x{Vhf (23) 

is the Hamiltonian of the Gaussian model. By identifying operators in the 
loop model with operators in the Gaussian model, the asymptotic behavior of 
correlations can be inferred by computing appropriate Gaussian averages. 

Unfortunately, there is no simple representation of (f> (L, z) in terms of the 
field h(r). However, there exists a related operator, defined on the infinite 
plane, that can be represented as a two point function. For two points X\ 
and x 2 in the plane, consider the operator <j) (z, x±, x 2 ) that gives loops winding 
round x\ or x 2 , but not both, a weight of z. All other loops have a weight of 1. 
Interpreting (f> (L, z) as a one point function, the exponents associated with the 
large L behavior of (<fi(L,z)) should be half those of ((f) (z, x\, x 2 )) describing 
the large |a;i — x 2 \ behavior. The operator </> (z, x%, x 2 ) can be expressed as: 

(/) (z, x ll x 2 ) = exp (i (0 + tt/3) h (xi)) x exp (-i (9 - tt/3) h (x 2 )) (24) 

Here 9 = arccos (z/2). An oriented loop winding round x± or x 2 , but not both 
points, is given a phase factor of exp (±i#) by including this phase factor and by 
dividing by exp (±i7r/3), which is the phase factor for any loop in the loop model. 
Loops winding round both X\ and x 2 pick up a phase factor of exp (±2i7r/3). 
This has the effect of changing the sign of the phase of such loops and thus 
leaves the weight unchanged. 

It follows from |13| that the exponents r„ associated with <j)(z 1 xi,x 2 ) are 
given by: 



3 [9 + 27m) 2 1 



(25) 



We thus expect that: 



>(L,z))= V A n (L,9)L-i^^+* (26) 



n=—oa 



Where the coefficients A n (L, 9) are analytical corrections and have the expan- 
sion: 

A n (L,9)=a n , Q (0) + Y / ^fP- (27) 



k=2 
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Table 4: Exact values of the amplitude function ao,o (9). 



9 


ao,o (#) 




even L 


odd L 





2i(!)*exp(-C'(-l)) 


I^2*(|)^exp(-C'(-1)) 


tt/3 


1 


1 


2vr/3 


2*(l) A ^t|exp(-C'(-l)) 


^(D^^expK'C-l)) 


7T 




? 


4tt/3 









1.2 




Figure 4: The function ao,o (#)■ 



The leading asymptotic behavior is thus: 

{(j>{L,z)) ~a ofi (9)L~l£ + ^ (28) 

Assuming that the a n ^ (&) are analytical functions, it follows from l|26|) that 
these functions for fixed k are analytical continuations of each other: 

a n ,k (0) = a , k (9 + 27m) (29) 

It also follows from l|26|l that: 

a_„, fc (9) = a n:k [-6) (30) 

By calculating the asymptotic expansions of the exact evaluations of I? (L, 9) 
for 9 a multiple of 7r/3 given in Table |21 it is possible to obtain exact values 
for function ao,o (9) at those values of 9 for both even and odd L. In appendix 
|I|we present a simple method to obtain the expansions. In Tabled we have 
listed the exact values for function ao,o (#)• We have also numerically computed 
this function for many values in the interval < 9 < lix. A graph of ao^ (9) 
is shown in figure 01 We can see from Table 01 that for odd L the amplitude 
for 9 a multiple of 7t/3 except tt is 2 C os(l /2) times the amplitude for even L. 
Numerical computations appear to indicate that this relation holds for all 9. If 
this is indeed the case then it follows from i|ltjfl and l|54|) that the determinant 
D (L, 9) has the same asymptotic behavior for even and odd L. 
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The function D (L,ir) is zero for odd L. The limit 



^)=lim D ^ a 

«-»*2 cos (6>/2) A HT (i) 

is, however, nonzero and finite. The divergence of the amplitude function at 
6 = 7r suggests the presence of the logarithmic term log (L) in the large L 
behavior of <fi (L). According to (|26|l all exponents should be 2/3 plus an integer. 
We thus expect that the large L behavior is given by an expansion of the form: 

4> (L) = (f; a k L-i- k j + (jr b k L-?- k ) log (L) (31) 

\fc=0 / \fe=0 / 

where a\ = b\ = 0. We have numerically verified this expansion. 

5 Enumerations of families of nonintersecting 
lattice paths and cyclically symmetric plane 
partitions 

In it is shown that the determinant 



det 

Kr,s<L 



r + s-2 
r-1 



exp (i9) S r! , 



(32) 



in H17J) gives a weighted enumeration of certain types of nonintersecting lattice 
paths and cyclically symmetrical plane partitions. Consider a family of nonin- 
tersecting lattice paths on an L x L square lattice subjected to the constraints: 

1. A path is allowed to move one step to the right or downward. 

2. A path is required to start from a point (0, k) with < k < L — 1 

3. A path starting from (0, k) is required to end at the point (fc, 0) 

See Fig.[S]for an example. The determinant l-'il'l) gives a weighted enumeration 
of all possible lattice paths satisfying the above constraints, where a configura- 
tion with s paths is given a weight of exp (i9s). 

Families of lattice paths of this type are in bijection with cyclically symmet- 
rical plane partitions [111 IIP) . A plane partition of an integer N is an array of 
integers n jjk , such that n j+1 , k > n jjk , nj^ k+1 > n hk and N = J^JLi J2T=i n 3,k- 
Plane partitions can be represented by a pile of unit cubes by introducing x,y,z 
coordinates in Z 3 and placing at position (J, k,0) a stack of n^ k unit cubes. A 
cyclically symmetric plane partition is a plane partition whose representation as 
a pile of cubes is symmetric under a cyclic permutation of the x,y,z coordinates. 
The bijection maps lattice paths of the above type to cyclically symmetric plane 
partitions that fit in an L x L x L box . Families of lattice paths containing n 
paths are mapped to cyclically symmetric plane partitions with n cubes on the 
main diagonal. 

Numerical calculations indicate that the expansion (|26(l is also valid for com- 
plex 0. It is thus possible to obtain the asymptotics for real values for the weight 
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E 4 
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Figure 5: A typical configuration of nonintersecting lattice paths enumerated 
by the determinant in (|32|l for L = 4. Paths starting at Ai have to end at Ei. 



of a lattice path from this expansion. Assigning a weight of w to a path yields 
for the asymptotics of the weighted enumeration TV (w) : 



Another interesting problem is the enumeration of nonintersecting lattice 
paths with fixed numbers of paths. The number of families with n paths is the 
coefficient of exp (in9) of the determinant l|32|l . which is the coefficient C n (L) of 
the characteristic polynomial of the LxL Pascal matrix, see (@J . It is possible to 
derive an asymptotic formula for this number if n = L/2. Fourier transforming 
the asymptotic formula for D (L, 9) w.r.t. 6 by expanding ao,o (&) about 9 = 
yields: 



Here we have used the exact value for ao,o (0), given in Tabled 

6 Probability that a point is on a non-contractible 
loop 

For even L loops can wind round the cylinder, while for odd L there is a path 
spanning the length of the cylinder. In |19| we have conjectured that the prob- 




(33) 




(log 



(34) 
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ability, W (L), that a point is on such a non-contractible loop is given as: 



wm =£w (35) 

In this section we will show that this follows from the conjectures (flip and 
1)14(1, Consider the two faces of the lattice next to some point on a loop. The 
loop can surround either one of the faces or for even L wind round the cylinder, 
or for odd L, it can span the entire cylinder. The probability W (L) is thus the 
probability that the loop does not surround either one of the two faces. Denote 
the probability that the face to the left is surrounded by n loops and the face 
to the right by m loops as P (L, n, m). Since both faces are surrounded by zero 
loops if and only if the loop does not surround any of the faces we have: 

W (L) = P (L, 0, 0) (36) 

Either the number of loops surrounding the left face differs by 1 from the number 
of loops surrounding the right face, or both faces are surrounded by zero loops. 
It thus follows that the probability that a single face is surrounded by m loops, 
P (L, m), can be expressed in terms of the P (L, a, b) as 

P (L, m)=P (L, to, to + 1) + P (L, to, to - 1) (37) 

for 1 < to < [L/2\, 

P(L,0) = P{L, 0,1)+ P(L, 0,0) (38) 

and 

P (L, [L/2\ ) = P(L, [L/2\ , [L/2\ - 1) (39) 

These equations, together with P (L,rt,m) = P(L,m,n), which follows from 
the fact that the loop model is symmetric under half turn rotations, imply: 

L/2 

W{L) = Y J P{L,k){-l) k (40) 

fe=0 

According to (|15l) this is (cf> (L, —1)) and according to conjectures (|llf> and 114|) 
this is given by D l L '' 2 ,Z(f > which simplifies to . A ^}^ (see Table El- 

7 Some conjectures about probabilities involv- 
ing connectivities of points on a row 

In this section we give a conjecture for the probability that n consecutive points 
on a row are on different loops and, for n < 4, conjectures for the probability 
that they are on the same loop. These conjectures are valid for even L. 

Let S (L, n) be the probability that n consecutive points are on different 
loops. We conjecture that: 
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where 



T(L,n) = 

mJ [L 2 -(2k + iy 

The limit L — > oo yields for large n: 



n-l-2fc 



(42) 



" s ' ( x ,;) x I 27 



(43) 



Let R(L,n) be the probability that n consecutive points on a row are on 
the same loop. We conjecture that: 



R(L,2) = 
R{L,3) = 
R(L,4) 



111 2 



16 (L 2 - 1) 
53L 4 + 44L 2 + 128 
128 (L 2 - l) 2 

5 1993L 8 - 11544L 6 - 12528L 4 - 150784L 2 - 294912 
'215 (i 2 - 9) (L 2 - l) 3 



(44) 



8 Conclusion 



We have presented new conjectures on exact expressions for correlations in the 
dense 0(1) loop model on L x oo square lattices with periodic boundary con- 
ditions. We have generalized an earlier obtained result for the probability that 
a point is surrounded by m loops to odd L. Using this conjecture we have cal- 
culated the expectation value of the operator that changes the weight of loops 
winding round a given point. It turns out that the expectation value is given by 
a binomial determinant that is known to give a weighted counting of cyclically 
symmetric plane partitions and also of certain types of families of nonintersect- 
ing lattice paths. By calculating the asymptotics of the expectation value of this 
operator, we have obtained the asymptotics of this binomial determinant. For 
special values of the loop-weights the determinant can be calculated in closed 
form. For those values we have exactly evaluated the amplitude function of the 
leading asymptotic behavior of the expectation value of the operator. 

Finally we have given conjectures on the probability that n consecutive 
points are all on different loops and, for n < 4, that they are on the same 
loop. 
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A Asymptotics of correlations 



The asymptotic behavior of correlations that are given as simple product formu- 
lae can be obtained by using the Euler-Maclaurin summation formula. However, 
this approach will yield the result up to a constant of proportionality. To obtain 
the results presented in Table El we need to find the asymptotics of A (L) and 
Aht (L), including the constant term. Both A(L) and Aht (L) are given as 
a product over factorials, and for such cases a simple method exists to obtain 
the asymptotic behavior, which is particularly suitable to be implemented in 
computer algebra systems. 

Products over factorials can be expressed in terms of the Barnes G-function 
0] 03 El EJ This function satisfies the relations: 

G(z + l)=T(z)G(z) 

G(l) = l (45) 

Where T (z) is the gamma function. The Barnes G-function is uniquely defined 
if the condition of convexity |2fi| is added: 

d 3 

log (G (z))>0 for z > (46) 

az J 

In the following we shall write xl for r (as + 1), also for noninteger x. It follows 
from the definition of the G-function that: 

n/i m G (n + 2 + x) , „„, 

( fc + *) != G (l + * (4?) 



fc=0 



The products occurring in A (L) and Aht (L) can be rewritten as in the above 
equation by using Gauss' multiplication formula for the gamma function. E.g. 
the numerator in A (L), see ©, can be rewritten as: 

3j+i j j j'-i 

(3j + 1)! = n fc = n 3fc n ( 3fc + v n ^ + ^ ^ 

k=l k=l k=0 k=0 

By taking out factors of 3 in the products this becomes: 

From (|4Tj) it thus follows that: 

TTro/HV- 3f " M " G(n + l)G(n+i)G(n + §) 

M ( J ) ' "(-D ! "K)<" Gf(f)G(|) (50) 

To obtain the asymptotic behavior of our correlations, we thus need the 
asymptotics of the Barnes G-function and some special values of this function. 
According to pQ, the asymptotic behavior of the Barnes G-function is: 

log (G (2 + 1)) =z 2 Q log* - |) + i log (2tt) « - 1 log (z) + C' (-1) 

" (51) 

\ - i>2fc+2 



^ 4fc(fc + l)z 2fe 
fe=i v ; 
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Here the B^ are the Bernoulli numbers, and Cf (z) is the derivative of the zeta 
function. 

The denominator of (|50"|) contains the term G (4/3) G (2/3). Such products 
can be simplified using the multiplication formula for the G-function |22J: 



G (nz) = exp [(l - n 2 ) (' (-1)] n « 2 - 2 /2— +5/12 

x{2n) ^-n*)Y[l[ G ( z+ i±l\ (52) 
i=0 j=0 ^ ' 



To express the product G (4/3) G (2/3) in terms of elementary functions one 
chooses z = 1/3 and n = 3. 

Proceeding in this way we have obtained: 



And 



(54) 



where 



r(i)/ -V3 



A R t (L) ~ H (L) for even L 
A HT (L) ~ 3-iH (L) for odd L 
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